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We investigate the asymptotic behavior of the Selberg-hke integral 




as — 7- oo for different scalings of the parameters a and b with A^. Integrals of this 
type arise in the random matrix theory of electronic scattering in chaotic cavities sup- 
porting channels in the two attached leads. Making use of Newton's interpolation 
formula, we show that an asymptotic limit exists and we compute it explicitly. 
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I. INTRODUCTION 

In his famous 1944 paper-, Atle Selberg introduced and computed the integral 

JV 

SN{a,b,f3):= / n \xi-x,\^^l[xr\l-Xif-'dx, 

>^[0,l]^ l<i<j<N ^=l 

_ fi' r(a + j(3)T{b + j/3)r(l + (j + 
11 r(a + 6+(iV + j-l)/3)r(l + /3) 

in the aim to solve a problem of Gelfond^. Since the sixties, many generalizations and 
applications have been developed (see^ for an interesting review about Selberg integral). 
The scope of these investigations involves many areas of mathematics : random matrices^"^, 
calculations of constant terms (see e.gr>^)^ symmetric functions^'^ (in particular. Jack and 
Macdonald polynomials^^ii), multivariate orthogonal polynomials^^, the value distribution 
of the Riemann C, function on the critical line^ among other applications. Since the seventies, 
many applications in physics have been found, especially in the theory of quantum Hall 
effeclr^'i^. The number of variables involved in the integrand is then interpreted as the 
number of particles. In this context it is interesting to study what happens when this 
number becomes very large. 

In the field of random matrices, the integrand of Selberg's integral corresponds also 
(for quantized values of /3) to the joint probability density of eigenvalues of one of the 
classical random matrix ensembles, the Jacobi ensemble. Matrices from this ensemble can 
be generated in three different ways: 

• as truncations of Haar orthogonal, unitary or symplectic matrices^-. For the case of 
unitary matrices, an important application arises in the theory of electronic transport 
in mesoscopic systems at low temperatures^, as detailed below. 

• as composition of Wishart matrices, with applications to multivariate statistics^. 

• as composition of projection matrices^. 

In the theory of quantum transport through mesoscopic devices (Landauer-Biittiker scatter- 
ing approachi2.i22i2i) ^ the wave function coefficients of the incoming and outgoing electrons 
in a cavity are related through the unitary scattering matrix 5* (2A^ x 2N , if N is the number 
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of electronic channels that each lead supports): 




(2) 



where the transmission (t, t') and reflection (r, r') blocks are (A^ x A^) matrices encoding 
the transmission and reflection coefficients among different channels. Many quantities of 
interest for the experiments are linear statistics on the eigenvalues of the hermitian matrix 
tt\ i.e., are quantities of the form A = Xlili /(^j) where /(x) is a smooth function (not 
necessarily linear) and Tj are the eigenvalues of tt^ (real numbers between and 1), which 
have the intuitive interpretation as the probability that an electron gets transmitted through 
the i-th channel. For example, the dimensionless conductance and the shot noise are given 
respectively by G = Tr(tt^)^° and P = Tr[tt''"(l — tt''')]^'^. The random scattering theory 
models the scattering matrix S for the case of chaotic dynamics as a random unitary ma- 
trix uniformly distributed in the unitary group, i.e., it belongs to one of Dyson's Circular 
Ensembles. 

From this information, the joint probability density of the transmission eigenvalues {Tj} 
of the matrix tt'^ , from which the statistics of interesting experimental quantities could be in 
principle derived, is readily recognized as the Selberg integrand (Jacobi measure) with 6=1 
and 2/3 = 1,2,4^- depending on physical symmetries of the Hamiltonian (for recent results 
on the use of Selberg integral in the quantum transport problem, see^-"^*). The importance 
of linear statistics, and their asymptotical properties when the number of channels grows to 
infinity, provides one of our main motivations for the present study. 

This paper is the continuation o^ and we are interested in the following integrals 



denoting averages over the Jacobi probability density. Especially when / is a power sum 
(section HTBD 



or a Schur function (section III AD . we raise the question of its asymptotic behavior for 
N ^oo. 

More precisely, using classical identities on symmetric functions we show that this problem 
reduces to the calculation of an inverse binomial transform. 





i<j i 
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This paper is the first step towards a combinatorial interpretation of the asymptotic behavior 
of Selberg-like integrals. Note that J.-Y. Thibon with one of the authors have already 
investigated the links between combinatorics and Selberg integrals^^— . 

The plan of the paper is as follows. In section Ull we give an expression of the integrals 
as a rational function in the numbers of variables A^. In section llllt we investigate some 
properties of the binomial transform which will be used to compute the limit values in section 
IIV[ Section |V] deals with several special cases related to combinatorics. Finally, in section 
IVII we provide concluding remarks, and in appendix |X] an alternative approach to corollary 
[ILl 



II. SOME SELBERG-LIKE INTEGRALS 



A. Schur functions 

Here we give an expression of ^^^^^^/'r^^'*^"''' , where 

sxix,, ...,Xn) = -pf (r'-r) 

L ^l<i<j<N\•^^ -^JJ 

denotes a Schur function (see e.g.-^). 

A rather classical formula which can be found in the book of Macdonald^(this is a special 
case of the exercice 7 p385) gives 



^s,^x„ . . . , x^)U = [[{X. -X,+J-^)l[ riK + a + b + 2N-^-l) 

i<j i=l ' 



Hence, 



(Sa(Xi, . . . ,XAr))a,fe _ j-j- Aj - Xj + J ~ ^ 



-A- r(Ai + a + - i) r(a + 6 + 2A^-i- 11 
X M 



f=T r(a + A^-i) r(Ai + a + 6 + 2A^-z-l)' 
First, we remark that the number of factors of 
1^ r(Ai + a + - i) r(a + 6 + 2A^ - i - 1) _ TT TT a + iV - z + j 

11 Tin J- N - A fT\~T17Tl^lI'9lv^^7^^r\ ~ 11 11 



J- r(a + A^-i) r(Ai + a + 6 + 2A^-i-l) i-i- i-i- a + b + 2N - i + j - 1 
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depends only on the partition A and not on the number of variables A^. In the same way, 
one has 



yr Xj- Xj+ ] - I -i-r Xi+ ] - I 
11 j — i j — i 



i<j i=l 

One needs the following lemma which is obtained by reorganizing the factors and simplifying 
the resulting expression. 

Lemma II. 1 For N large enough^, one has 

N 



nb + i _T-ra + 6 + i 
b + c + i ~ ^^b + N + i 

=a+l 1=1 



where a, c G N. 



Proof First write 



n 



b + i 



n 



b+i+a+l 



b + c + i J--!- b + c + i + a+1 

c-1 N 



Y[ib + i + a + l) Y[{b + i + a+l) 



N-a-c-l 



N 



n 



i=0 



But 



and 



TV 



■ c + I + a - 



N-c 



1) Yl (b + c + i + a+l) 



i—N — c—a 



Y[{b + i + a + l) = Y[{b + c + i + a + l), 



i=0 



TV 



c-1 



Yl {b + c + i + a + l) = Y[{b + N + i + l). 

i=N—c—a i=0 

Hence, by substituting these two identities in ([3]), one recovers the result. □. 
If one applies lemma ULT] to njL^(A)+i ^^j^T^' finds 

TV . . A, . 

-Q Aj+j-z _-Q J + N -I 



1 — I 

j=£(A)+l j 



Hence, 



(3) 



^ The condition 'W large enough!^ can be omitted provided that we use the notation ni=a /(* 
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Proposition II. 2 One has 

{sx{Xi,...,XN))a.b 



n 



n -"^W^xi 



j=i+i 



3 - * 



{j + N -i + l){a + N -i + j) 
(^(A) +j -i + l){a + h + 2N -i + j -I) 



If a = aiN + ao and 6 = biN + 6o ore two linear functions of N , this implies that 



{Sxjxi, . . . ,XN))a,b 

where \X\ = A^. 



1 + ai 
Af^oo V 2 + ai + 6i 



|A| eW 

n 

1=1 



j=i+i 



J=0 



(4) 



B. Selberg-like integrals with a power sum in the integrand 



In this section, we study the integral 



h :-- 



{Pkjxi, . . .,XN))a,l 



where pk = Y^iX^. 

One uses the formula (see e.gM-'^): 

k~l 

Pk = y^(~i)'-s[(fc-»)i']) 

i=0 

where [{k — i)V] denotes the partition [{k — i), l, ■ -j ^^j - Hence, 

{^[{k-i)V])a,b 



k-l 



i=0 



From proposition 111.21 one has 



Corollary II. 3 For each k > 0, one has 



J 1 .y(k-lV-^ \N + j){a + N + j-l) 

^ h\ 2^^ M i j 11 a + b + 2N + ] -2 ' 



i=0 



j=-t 



Proof From Eq (j4]), one has 



{Slik-^)l^])a,b _ '-^ {p + N){a + N -1+p) 



n^^-^ n 



p=2 



^Li {i+p+l)(a + b + 2N+p-2) 



n 

J=2 



first part of tlie partition 
{N-j + l)ia + N-j) 
ii^j + 2){a + b + 2N-j-l) 



other parts of the partition 



(5) 



(6) 



(7) 
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But, 

n*"*""^ k — i + p — 2 V-i- k — i + p f k — 1 
p-l "J-l p+i ~l i 

p=2 ^ p=0 ^ ^ 

n (z+p+i) n (^_^-+2) " n ~n~=;fc!' 

p=0 ' ^ ' ^ j=2 ^ ' p=i+l ^ j=\ 

and 

y^ (iV-j + l)(a + iV-j) -pi (iV + j)(a + iV + j-l) 
Using these equalities in ([7]), one finds 



p.-i)i»])a,6 1 /A^-lV-pT-' (Ar + j)(a + iV + j-l) 



(l)a,6 A;! V ^ y n. a + 6 + 2iV + j-2 

The result is obtained by replacing ^ [(fe-^)^'^^"''' with its value in (I6]).n 

The expression in corollary III.3I is the starting point in computing the limit lim7v->oo 77 
(in section [IV|) . using the tools we are going to introduce in the following section. 



III. INVERSE BINOMIAL TRANSFORM 

A. Inverse binomial transform and Newton's interpolation formula 

In this section, we shall use widely the inverse binomial transform (see e.gM^) that 
operates on a sequence of polynomials F = {fi{x))i by 

Proposition III.l Let F{y) = Yl^=o '^ii^)y^ unique polynomial in y with coefficients 

in C[x] of degree p (in y) interpolating the points 

(0,/o(x))...(p,/p(x)). 

We have: 

p 

^,'[¥]=k\J2Si,kai{x), 

i=k 

where Si^k is a Stirling number of the second kind. 



7 



Selberg-like integrals 



Proof By linearity, it suffices to show the result for F{y) — yP. We remark that in this case 

i=Q ^ ^ 

by means of the well known formula 

i=0 ^ 



□ 

Consider the divided difference operator dy^^y^ acting on the right of any expression / in 
{yi,y2} by 

J'^yiy2 — J 

where (Jy^y^ permutes yi and y2 in f. 

Reiiicirk III. 2 Note that, assuming that k < p, the degree of yody^y^ . . . dy^_^yi^ is equal to 
p — k. Therefore, ifg{x,y) is a polynomial of degreep inx andy, the degree oftB^^ [{g{x,i))i] 
equals p — k. 

This operator is the main tool to describe the Newton interpolation. Indeed, consider a 
one variable function f{y) and a set of interpolating variables {yo, . . . , y^}- One has 

f{y) = f{yo) + f{yo)dyoyi{y - i/o) + • • • 

• • • + f{yo)dyoyi ■ ■ ■ dy^_,y^{y -yo)---{y- yk-i) + R{y) 

with R{yi) = for each i = 0, . . . , k. 
We will denote 

fdm...n — f{ym)dy^y^_^_j^ ■ ■ ■ , ^y^y^_^_^\y^=i, 

for each pair of integers m < n, with the special case fdm...m = /(^)- Set, also, di := 
With this notation, the polynomial / of degree n — m interpolating the points 

(m,/(m)),...,(n,/(n)) 

becomes 



f{y) = Y.f^rn..,{y-m)...{y-{j-i)). 



j=m 
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Remark that the Stirhng numbers appear when we write in terms of falhng factorials, 
{y)k ■.= y^y-l)...{y-k + l), 

p 

yP = ^Sp^k{y)k- 
k=0 

This means that the Stirling numbers are the coefficients in the Newton interpolation of 
at y = 0, . . . ,p. With our notations this reads 

Sp,k = y''do,„k+i. (9) 
Hence, by linearity, one obtains immediately the following result from proposition IIII.ll : 

Corollary III. 3 Let k < p be two integers. Let F be the unique polynomial in y with 
coefficients in C[x] of degree p (in y) interpolating the points 

(0,/o(x))...(p,/p(x)). 

We have: 

03^^ [F] = k\Fdo,„j,.D 
When acting by do___k+i on y^, one observes the following (shifted) induction. 



Proposition III. 4 



y^do...k+i = iy + ir-% 



Proof Since 



p-i 

yo^yoyi ■ ■ ■ ^VkVk+i ~ yi ^yiy2 ■ ■ ■ ^ykVk+iyoy 
1=0 



the specialization gives 



y^do...k+i = y^ ^di...k+i. 



We conclude the proof by noting that the coefficients in the Newton interpolation of any f{y) 
aty = l,...,/c + l are, respectively, equal to the coefficients in the Newton interpolation of 
fiy + 1) aty = 0,...,kn 
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B. An example of generalized (inverse) binomial transform 

As an application, consider the polynomials 

k—i—l i—1 

i=o i=o 

For simplicity, we will denote := Pi{x] a, h) when there is no ambiguity. 
Proposition III. 5 When p < k, one has 

p— 1 k—p—l 

<B,'_^[P^\...,P^] = ll{x + b-t) n (b-a-p-t) 

i=0 i=0 

Proof First remark that Z/o^yoyi • • • is a symmetric polynomial in {?/o, • • • , 

Hence, we can permute the variables in the expression and obtain 



yovi 



Applying the same argument to i/jdy.y^ -^ . . . dy^^yg which is symmetric in the variables 
{yo, . . .,yj} one gets 

^i'^ww-i ■ ■ ■ ^yiyo^yoyj+i ~ Vo'^yoyi ■ ■ ■ ^j/j-i%^s/o%+i- 
By definition of dy^^y.^^ one obtains 

ip. a _ yj+i'^yj+iyi'^yiyi ■ ■ ■ ^Vj-^Vj ~ Uo^yoyi ■ ■ ■ ^yj-iyj /-.^.n 

Vj+i - yo 

Again, 

yj+i^yj+iyi^yiy2 ■ ■ ■ ^Vj-iVj — Vi^yiyi ■ ■ ■ dy^yj^^, 

and eq. (ITUI) becomes 

ip) rj _ yi9yiy2 ■ ■ ■ 9y^yj+i ~ Vo^yoyi ■ ■ ■ ^yj-Wj (-,-,\ 

yo'^yoyi ■ ■ ■ '^i/ji/j+i ~ 0,0, ' ^ ^ 

Vj+i — Vo 

Let us prove the result by induction on k. Note that ii k = p, the result is straightforward. 
Denote by P{y) the unique polynomial of degree k — p + 1 in y such that P{i) = Pp_^^ for 
each i = . . . k — p. By linearity eq. (fTTj) gives 

p/ \ _ P(yi)^t/lt/2 ■ ■ ■ (^Vk-p-lVk-p ~ ^{yn)dyoy;^ ■ ■ ■ t^t/fc_p_2^fc_p_i 

yk—p yo 
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Specializing at i/i = i, one obtains 



Pf^i...fc-p ~ P<9o...fc-p-i 



Pdo...k-, = "-""^ (12) 
k — p 

By definition of P and using corollary IIII.3[ one has 

"^k-piPp, ■■■,Pk] = ^-'[Pio), ...,p{k-p)] = {k- py.Pdo...,.,. 

Hence, eq. (fT2i) yields 

^,'_^[P,\ ...,Pl:] = {k-p- i)!(Pai...fc_p - Pdo...k-p^i). 

And then, 

^fe-pl-^p ' • • • ' ~ ^fe-p-il-^p+i' • • • ' ~ ®fc-p-i[-^p 5 • • • ' -Pfc-i]- (13) 
Remarking that ii i < k one has 

P''=(a;; a, b) = {x + a)P^-\x- a + 1, 6) 

and if i > 

a, 6) = (a; + h)P^-\x- a, 6 - 1), 

eq. (|T3l) becomes 

»,7!p[Pp^ . . . , P,fc] = (x + [^P a, 6 - 1), ... , P;fri'(.T; a, b - 1)] 

One recovers the result by applying the induction hypothesis. □ 
Note that, setting p = in the previous proposition, one obtains 

k-l 



^-,'[Po^...,Pk'] = ll(^-a + ^). (14) 



1=0 



Remark III. 6 Remark \III. B and proposition \III. 5l imply that the (Newton) polynomial in- 
terpolating the point (0, Pq), . . . ,{k, P^) (i.e. the polynomial P of degree k in y such that 
P{x,y) = P^) is a bi-variate polynomial (in x and y) whose degree is k. 
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C. Leading coefRcients 

Let i{x,y) be a bivariate polynomial of degree m. In this section, we investigate the 
coefficient of the leading term of the binomial transform 53^^[(f(x, 2)^^)4], that is the scalar 



yi=i' 



multiplied by k\, where [x*]P(x) = if P{x) = Yli'^i^i- The following result explains how 
to manage the shift induced by the multiplication by i'^: 



Proposition IIL7 One has 



otherwise, 



where fp(x, y) = f{x,y + p) . 

Proof Proposition IIII.4I yields 

f{x,y)y^do...k+i = f{x,y + l){y+lY-%...k. 

Note that, since f (x, y){y + 1)^ is a bivariate polynomial of degree p + m — 1, 

[x^^"'-']f{x,y + l){y + ly-' = [xP+^-'']f{x,y + l)^/)^-^ + R{y), 

where i? is a polynomial of degree at most k — 1. Hence, since the operator dy^^y-^ ■ ■ ■ dy^_^y^, 
lowers the degree of k in the ?/j's, one observes 

[x^+'^-^]i{x,y + l){y + ir~^do...u = [x^^^^-^^x^y + l){yy-^do...k, 
or, equivalently, 



£,_i,p_i(fi) ifl<A; 
otherwise, 



Iterating the process, one proves the claim. □ 

If {ai{x))i is a sequence of polynomials, one defines 
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Remark III. 8 From remark [III. (H if i(x.y) is a bivariate polynomial of degree m, its trans- 
form, 1l'^[i{x, 0), . . . , f (x, k)], is a polynomial in x whose degree is also (at most) m. Indeed, 
the polynomial F of degree k in y interpolating the points (0, /(x, 0)Po ), . . . , (fc, /(x, k)P^) 
is a bivariate polynomial whose degrees in x and in {x,y} are m + k. Acting by 9o...fc, one 
obtains a polynomial in x which is a linear combination of coefficients [y^]F (whose degree 
in X is m + k — i). We conclude by noting that, when i < k, [y^]F gives no contribution since 
Vo^yoyi ■ ■ ■ ^Vk-iVk '^^ ^ polynomial of degree i — k. 



In particular, when f(x, y) = y^, eq. f lT^ and proposition IIII.7I imply 



[x^]^fW) 



k\ 



k—p—l 



{k-p)\ 



Y\ {b — a — p — i) ii p < k 

otherwise. 



1=0 



(15) 



The following result allows us to compute the leading term in the action of T^''' on a product 
of linear factors of the form ex — dy + e. 



Corollary III.9 



a,h 



fp-i 



k I 



E 



k-j-l 



\ {o,...,p-i}= I 1=0 

\ {si,...,Sj}U{ti,...,tp_j} 



Proof First note that the variables give no contribution: 



[xP]%k 



Y\ {cjX — dji + Cj] 

Kj=0 



[xm. 



Y\ {cjX — dji) 
\i=o 



Indeed, the coefficient of Cj^ . . . cj^ in 11^=0 i^j-^ ~ ^iV + ^i) ^ polynomial /(x, y) whose 
degree is p — s. From remark [Tll.2t the degree of [(/(x, i))^ being at most p — s, it follows 
immediately that [x^'j'Tfe [{f{x,i))^ ^ only if s = 0. 
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Now, to obtain our result, it suffices to expand YYj^l i'^j^ ~ d-jv) as a polynomial in x and 



[^1 E(-i)''' 

j=0 

j=0 



p-1 

n (^^'^ ~ + 

j=0 

I 



51 • • • Cgjdt^ . . . dtp_ 



( 



1 {0,...,p-l} = 
\ {"1 Sj}U{ti 'p-j} 



. . . Cg.dt^ . . . dt^ 



1 {0,...,p-l} = 

\{si,...,Sj}U{ti,...,tp_j} 



Applying eq. f|T5|) . one recovers our result. □ 

In the next section, we use the results of this section in the aim to investigate the asymp- 
totic behavior of the integrals ^. In particular, we show that its convergence (section [IV]) 
is a direct consequence of remark IIII.8I and we compute explicitly the limit (section IIVBII 
by means of corollary IIII. 91 



IV. ASYMPTOTIC BEHAVIOR OF k 



N 



In this section, we use the tools described in the previous section to prove the convergence 



of the integral ^ and compute the limit. 



N 



A. Convergence 



Suppose now that a = a{N) and b = b{N) are linear function of A^. One has: 
Theorem IV. 1 

4 

lim — < +00 

N^oo N 

Proof We start from Corollary III. 31 and write 

4 1 ^fc(A^) 



N klNYl'Zlk+MN) + b{N) +2N + J-2) 



where 



k-l 



k - 1 



Y[ {2N + a{N) + b{N)+j~2)^, 



fe-i-i 



k-l 



n (N + j){N + a{N) + J - 1) n {2N + a{N) + b{N)+j-2). 

j=-i j=k-i 



(16) 
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We need the following lemma: 

Lemma IV. 2 The degree in N of the polynomial ^k{N) is 2k. 
Proof For convenience, set a = aiN + oq and b = biN + bo, and write 

-i-l k-1 

Y[ {2x + a{x) + b{x) +j -2) Y[ {2x + a{x) + b{x) + j - 2), = 

j=-k+l j=k-i 

k-i-2 i-l 

Yl ((2 + ai+ bi)x + j + ao + bo - I - k)Y[{{2 + ai + bi)x + ao + bo - j + k - 3). 

j=0 j=0 

With the notation of the previous section, one recognizes 

k-i-2 i-l 

Yl ((2 + ai + bi)x+j + ao + bo - 1 - k) ]J((2 + ai + bi)x + aq + 5o - j + fc - 3) = 

3=0 j=0 

Pt~^{{2 + ai + bi)x- ao + bo - k - l,ao + bo + k - 3). 
If one sets 

fc-i 

2 + ai + 6i ' / V 2 + ai + &1 



3=0 

the following holds 

Hence, from remark Ull.St the degree of ''XlkiN) equals the degree of Qfc(x,y) that is 2k. □. 
The degree in N of the denominator 

k-l 

N Yl {a{N) + b{N) + 2N + j - 2) 

j=-k+i 

of ^ is 2k. From lemma IIV.2t ^ is a rational fraction in N whose numerator and denomi- 
nator have the same degree 2k. Hence ^ converges. □ 

B. Computation of the limit 

Let F^^^^^.^^^^^(x,?/) = ITjZliaiX+j -y + (3i){a2X+j -y + (32). And set, for convenience, 
FP . (x) = (Fl . .„ . (x,i)) One has: 
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Proposition IV. 3 The coefficient of x^"^ in ''K^ ^^.q,2 ^2(2^)] does not depend on /3i and 
^2- More precisely, one has: 

Proof This equality is obtained from corollary IIII.9I setting cj = ai, cj+p = 02, dj = dj+p = 
1, Cj = j + (3i and e^+p = (32 for each j = 0, 1, . . . , p — 1. □ 
Using this result, one finds : 

Theorem IV. 4 Setting a = aiN + oq and b = biN + bo, one has 

Proof We have seen (see eq ( IT71) ) that 
lim — 



iv^ooiV A;!(ai + 6i + 2)2'=-i 

^ r T\T2klrr-a.o+bo-l-k,ao+bo+k-3 r-rw-k | 

2 + ai+bi' '2 + ai+bi' 

By proposition IIV.31 one obtains 

i=0 j=0 ^ > 

Ef k\ f k \ f l + fli ^ 

One recognizes 

And, slightly rearranging the terms of the sum, we show the theorem. □ 



jj \2k — i — j J \2 + ai + 6i 



Let us rewrite ffTSl) as 



lim ^ = 1+5^ ((i + a^)fe-i + 



. fc — i — 1 V V \2 + ai +bi J \i + j + I J \ j 



(19) 



For convenience, we set Oi = £1 — 1 and 61 = — £1 — 1. With this notation, equation ( IT9 
reads 
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Rearranging the sum, one obtains 

or equivalently, in terms of binomial transform one has 

V. SOME SPECIAL CASES RELATED TO COMBINATORICS 
A. Simplest cases 

Corollary V.l 1. If ai = —1 and —\ then limAr_j.oo 7^- = 0. 
2. If ai 1 and bi = —1 then limAr_j.oo 7^ = 1- 

Proof The assertion (1) is straightforward from (fT8l) . 

To prove the second assertion, we need the following lemma 

Lemma V.2 Let a, b and c be three integers and denote 
With this notation, one has 




Proof By induction, remarking that 

□. 

Under the specialization bi = —1 (or equivalently £2 = ^1) and using the notation of lemma 
rOl formula ([20]) reads 

But lemma EJ yields { = C'fc^) = 0- This implies our result. □. 
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B. Central binomial coefficients and the specialization ai = bi = 



Under this specialization, ( !T8|) reads 



But, one has 

^ ( _L T I 1 ) ( ) ~ ( I T I 1 ) ■ 



' ' ^ k \ fk\ f 2k 



1=0 

Indeed, formula fl22|) follows from a well known equality 

oo 

E 



well known equality 
a\ f b \ f a + b 



J / \c + J J \a + c 



which can be proved by a straightforward induction on b. Hence, from fl22l) . one obtains 



k-l 



M Okie 2-^ \ 9 
3=0 

We need the following lemma : 



lim , / 1 II 11, 

N^oo N 2^k \ 2 I \ j J \k + j + 1 



j + k-l\ f 2k 



Lemma V.3 One has 

n—m / , ■ \ / o \ I2n^ 



j J \n + m + j J \mJ (^^) 

Proof Set Fnj := (—2)^-' (^"^^^'' j (n+m+i)- compute the associated Gosper sequence: 

r - 2 + p 

m — n + j — 1 

and we check that the two sequences verify the equality 

(2n + l)F„j + (m - - l)F„+ij = Gnj+i - Gnj. 

It follows that 

, / n \ , , / 72 + 1 

2n + 1 ( ) = n - m + 1 ( 

\m/ \ m 

Hence, the result is obtained by induction. □ 



From rvT2l one obtains 
Corollary V.4 

lk^J_/k\^J_f2k\ 
N^oo N 2^k\l/ 2^>'\k )' 
This result is well-known in the quantum scattering setting in symmetric cavitiea^"— . 
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C. Catalan triangle and the specialization oi = 



Let us set for convenience bi 

fe-i 

in — = — — 

N- 



1. Under this specialization eq f[T8|) reads 

k-i-j 



j + k-1 
3 



1=0 



Using formula (1221) . we obtain 



lim 



fc-i 



+ 



A; 



2k 

k + j + 1 



or equivalently 



lim — = , 

N^oo N k{l+e) 



fc-1 



j=Q ^ 

1 k—l/k—l — i f k 1 

i+£)2fe-i 51 H ^^^y ( J 

1=0 V j=0 ^ 



+ j - 1 \ / 2fc 



By rearranging the factor of the products appearing in the coefficient of each one restates 
this expression in terms of inverse binomial transform: 

k-l 



{2k)\ 



-1 



^ ^ \k—l—i 

N^ooN k\{l + iy^-^ ^ l\{k - t - 1) 



We need the following lemma: 
Lemma V.5 

Proof First remark that 



1 



ij + k){j + l + k) 



(23) 



{p + i){p + i + 1) 



;-l)™(m + 1)! 



P-_,(p;0,m + l) 



Hence, 



53: 



{p + i){p + i + 1) 
We conclude by using eq. (IT^ . □ 



[Po'"(p; 0, m + 1) . . . P-(p; 0, m + 1)] 



nr=r(p+o 



Now using lemma IV.5I in equality ( 123|1 , one finds 
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Proposition V.6 



^ ^ k-i f2k 



k \ i 

lim 



The triangle D := (T))fcieN sometimes called Catalan triangle (see e.g. sequences 
A008315, A050166 and A039598 in^^). 

Note that these numbers are related to many combinatorial objects. For example, they 
appear in the expansion of odd power of x in terms of orthogonal Chebyshev polynomials 
Uk{x) of the second kind (see e.g.— p. 796), since 



k-l 

„2k~l _ 



i=0 



Another example is given by R.K. Guy in^: he showed that the number of walks in a lattice 
with k steps (each in direction N, S, E or W) starting at (0, 0) and at a distance i from the 
X-axis equals Dk+i-i,k+i- 

D. Symmetric Dyck paths counted by number of peaks and the specialization 

6i = 

For convenience, let us set Oi = i — This case has been already computed by M. Novaes 
in^i^. With our notation, he proved the following formula 

„,„ = + 1) f (zl)!l f * - 1 - m ' 

Our goal is to identify the coefficient ai^k such that 

lim Jl- 



N^ooN (l + £)2'=-l' 

Let us sketch the proof of the following result, announced by two of us wM. 



Proposition V.7 



h i '^^^ fk-l\ fk-1 

nAr-(l + £)2-i £^ Ull A [|J 
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Proof First rewrite eq as 

N^oo N 2k -I ^ i \ i-l ^ ^ > 

1=1 \ / \ / 

Expanding (1 + and rearranging the sum one obtains 



N^ooN (l + £)2fc-i ^ ' \ i-j J j + 1 

i — \j — 

Hence, we need to prove: 

=.^=D-..f'':/#=(,i,)(,;,)- 

Using Gosper algorithm and the Zeilberger method (see e.gM), we find that the sequence 
2k,i is completely determined by the following three terms relation: 

- {2k + l-i){2k-i) \i + {2 + 2i-2k) 2k,i+i + + + 2) 2k,i+2 = 0, (25) 

and the initial conditions 

\o = 1, = k. (26) 

k \ f k 



A straightforward calculation shows that the numbers p"| j y J verify also eq f l25|l and 
( |26|) . This concludes the proof. □ 

Remark that the numbers Jik^i have an interesting combinatorial interpretation since it 
is the number of Dyck paths of odd semi-length 2k — 1 with i peaks (see e.g.— and sequence 
A088855 in^S). 

VI. CONCLUSION 

This work is the first step towards a combinatorial interpretation for the asymptotic 
behavior of the integrals 

Wn-={xi---Xn)-= jrTc (27) 

for any partition A where 

{f{xi,...,XN))lb:= / f{xu...,XN)T\\Xi-Xj\'^''T\x'^-\l-Xi)'''^dXi. 
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In particular, our goal is to find an algebraic proof of the following conjecture suggested by 
numerical evidencej^. 

Conjecture VI. 1 One has: 

1 '^'^ 1 

i^L^iv^^"^"^^' • • • ' "^^^^^ = n i^iL^iv^"'^"^^' • • • ' ^^^^ 

i=l 

where p\ stands for the product px^ . . . p\^^^ ■ 

The limitations of the method described in this article is that we need an explicit expression 
for the integrals. This quickly becomes tedious in the general case. For testing equality fl28l) . 
we first obtain an expression of a product of power sums in terms of Jack polynomials, we 
implement an identity similar to (jl]) for the integral associated to each of Jack polynomials 
and finally we take the limit after simplifying the sum. 

The link between conjecture IVI.ll and the integral ( 1271) is the following. First, two of the 
authors remarked in^' that integral f l27|) can be restated in terms of monomial symmetric 
functions 



where P\{N) is an explicit polynomial whose degree is exactly ^(A). Remarking that 



mx = axPx + ^ a^Pf,, 



and assuming conjecture IVI.ll we see that lim7v^oo(A)5v equals (up to an explicit multiplica- 
tive coefficient) limAr^oo(]vi(ATPA(a;i, . . . , a; at))". 
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Appendix A: Calculation of (corollary III.3[) via spectral density 



We sketch here the derivation of an alternative formula for (Corollary III. 3P via spectral 
density of the Jacobi ensemble. 

The spectral density of the Jacobi ensemble is the marginal of the joint density function 
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(integrand of the Selberg integral), defined as: 

„1 JV 

Pn{xi; A,B) = N dx2--- dxN JJ l^i '^k^W^f'^- ^i)^ (Al) 

"^0 j<k i=X 

and it is normalized to A^, jj^ dxip]y{xi; a, b) = N. 

Such density is known analytically for all in terms of Jacobi polynomials P^^'^^ (^z) as: 

N-l 

Pn{x; a, B) = x^{l - x)^ B)[Pj'^'''\l - 2x)]^ (A2) 

j=0 



where: 

(2j + A + B + l)r(j + l)r(j + A + B + V 

r(j + A + i)r(j + 5 + i) 

The quantity Ik in Corollary III.3I is a linear statistics on the eigenvalues of the Jacobi 
ensemble. As such, it can be computed as a 1-fold integral over the density: 

dx x^ Pn{x]A,B) (A4) 



where ^4 = — 1, 5 = 6—1. 

The integral (lA4p can be evaluated in terms of nested finite sums as: 



h = Y. '-"^^ E dmM. B)f^^^{A, B)gtUA, B) (A5) 



j=o •' m,e=o 
where: 

{-j)m {A+B+j + l)m {A+m+l)j 



dmj{A,B) 

Ak)(A R\ ^ r{A+j+l)r{A+k+rn+l)r{B+l) 
JmjK^J-^) j\ r{A+l)r{A+B+k+m+2) 

Sk) ( A n\ _ (j+A+B+l)e {A+k+m+l)e 

ymjn^^-^J l\ {A+l)i {A+B+k+m+2)t 

where {x)n = T{x + n)/r{x) is a Pochhammer symbol. 



The identity ( lASp is obtained straightforwardly by first expanding one of the two Jacobi 
polynomials using the definition in^A and then computing the remaining integral using for- 
mula 7.392.1 in^-. The equivalence between formula (lASp and corollary III. 31 can then be 
proved by means of elementary but lengthy algebraic steps. 



25 



